First passage time distribution for a random walker on a random forcing energy landscape 
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We present an analytical approximation scheme for the first passage time distribution on a finite interval of a 
random walker on a random forcing energy landscape. The approximation scheme captures the behavior of the 
distribution over all timescales in the problem. The results are carefully checked against numerical simulations. 

I. INTRODUCTION 

The dynamics of a random walker on a one-dimensional random forcing (RF) energy landscapes has been a subject of much 
interest over the last three decades. Part of the interest is due to dynamics in quenched disordered systems. Examples are the 
dynamics of a random-field Ising model |Q-Q] and the motion of dislocations in disordered crystals Q]. More recently much 
interest has been due to many application in biophysical settings. It seems that in such systems random forcing energy land- 
scapes are the rule and not the exception. In this context, the dynamics of random walkers on random forcing energy landscapes 
have found applications in the mechanical unzipping of DNA Jfj-Q], translocation of biomolecules through nanopores |Q] and 
the dynamics of molecular motors Ilfj|4l3ll . In many of these experiments the first passage time (FPT) distribution is directly 
measurable and in some, such as in the translocation of biomolecules through nanopores, it is the most direct measurement. 

On a lattice the dynamics of the random walker is defined through the hopping rates between neighboring sites. We denote 
the hopping rate from site i to site i + 1 by pi and from site i to site i — 1 by g,. To generate a RF energy landscape the sets of 
rates {/?,} and {g,} are drawn randomly from distributions p (/»,) and q (qi) respectively (both distributions are independent of 
z). It is convenient to introduce an energy difference variables {£,■} such that £jp = e~ E ' (measuring energies in units of ksT). 
Thus each realization of the environment in the RF model can also be drawn as a set of i.i.d. random variables {g,} and {Ej}. 
Note that the energy landscape in this model is itself described by a biased random walk in energy space (see Fig. [T]for an 
illustration). 




d 



FIG. 1 : An illustration of a RF energy landscape. At site d + 1 there is a reflecting boundary. Here d — 50 and jl — 0. The energy difference 
between neighboring sites, Ej, is drawn from the Gaussian distribution. 
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It is well known IU4H17H that the long-time and infinite lattice asymptotics of a random walker on such an energy landscape 
is rather rich. In particular, the dynamics are controlled by a parameter, ji = \fi\, defined through the non-zero solution of the 
equation 

/ / n \ **\ 

= 1, CD 




where the angular brackets represent an average of realizations of disorder or, equivalently, over the distributions of {/?,} and 
{qi}. For jJ. > 2 the behavior is similar to that of a biased random walker on a flat energy landscape. Namely, for large times 
the mean position (3c) and its variance / (x — x) 2 \ both grow linearly in time (here the overline denotes an average over histories 
of the system starting from the same initial position). In this regime both quantities are self averaging. When < < 2 
the behavior is anomalous. For 1 < ji < 2 the mean displacement is self averaging and 3c ~ (3c) ~ t. The diffusion, however, 



behaves as (x — x) 2 ~ t 2 ^ while its average over realizations of disorder behaves as (jx — x)j ~ £ 3 ~^. For < < 1 both the 
drift and the diffusion are anomalous. Asymptotically the velocity, defined through lim (3c) /t, vanishes. Specifically, the drift 
behaves as x ~ (3c) ~ t 11 (note that although 3c and (3c) have the same scaling 3c is not a self-averaging quantity Jig! ). In addition 



(x — x) 2 ~ t 2 !^ and its average over the realizations of disorder behaves as \ {x — x)j <~ f 3 Finally, when /I = (commonly 
referred to as Sinai diffusion) 3c ~ In 2 1, (3c) =0,J?~ (x 1 ) -ln 4 f |Ql and (jc-3c) 2 ~ ((x-3c) 2 \ ~ t° 1 20]. 



In this paper we provide an approximation for the FPT distribution of a random walker on a RF energy landscapes in a finite 
interval. Namely, we are interested in the disorder averaged FPT probability density from site z'o > to the origin (site 0), 
F (t) — (f/ ->o ( f I {Pi j <?/})) with reflecting boundary conditions at lattice site d + 1 with d > iq, so that p e / = or = °° (see 
Fig. [TJ 1 . The results are summarized in Sec. [VT] Note that this rather rich behavior makes it impossible to write the averaged 
prop agator of the process as a scale invariant function, except in the Sinai's case fi = 0. Therefore, the techniques developed in 



[21 



2211 to calculate the distribution of FPTs for scale invariant processes are not directly applicable. Our results are compared 
with numerical calculations and shown to agree very well. To obtain the approximation we present a "random tilt" (RT) model, 
similar in spirit, to that used in 12311 . The parameters which define the RT model are given in terms of the parameters of a 
random walker on a RF energy landscape (for simplicity we assume a walker on a lattice) with no fitting parameters. This 
model naturally exhibits all of the regimes exhibited by other RF models after a proper disorder average. Specifically, within 
the model we obtain an exact result for the Laplace transform 

F(s)= f F(t)e- S 'dt (2) 
Jo 

of the FPT between two given points with specified boundary condition. The Laplace transform can be easily inverted numeri- 
cally. 

The paper is organized as follows: in Sec. |II]we present a brief review of relevant known results for a random walker on a RF 
energy landscape. In Sec. |III]we present and solve the RT model. In Sees. IV Bl and lV Al we determine the parameters of the RT 
model which are used to approximate the FPT probability density of the RF model. In Sec. IV Cl we compare the approximation 
to numer ical results. 

1 To obtain /o < with boundary condition cjd<i = results one simply takes 

i'o -> -k> 
H -> Pi 
Ei -> -Ei 

d -> -d. 
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II. REVIEW OF SOME KNOWN RESULTS FOR A RANDOM WALKER ON A RF ENERGY LANDSCAPE 

While the exact full FPT distribution, F (f ), for which we provide an approximation, is not known, several related results 
are known. Specifically, the mean FPT, defined through (t) = fntF(t)dt, can be obtained using the expression for a given 
realization of disorder (namely, a given set of {Ej} and {g, }) 12411 

io / i d-\ d-\ i k-i \ 

? =i ( 1 n-+i 1 n- > 0) 

where, as stated above, ^p- = e '. Therefore, the disorder average of t is 



4i/ {{e- E i)-\f \9i/ <«-*»>-! 



Note that this quantity may be very different from the typical FPT, defined through exp(lnf) (see for example [25] for a 
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discussion of the Sinai's case pi = 0). From Eq. (0 one may see that in the large d limit the leading order of (In?) is £ ln^ 1 = 
— Ziflf, where £ = (£,) is the average tilt of the RF energy landscape. Therefore, the typical FPT scales as 

exp(lnr) ~ e ~ Ed , (5) 

increasing exponentially with d for E < and remaining constant for E > 0. 

In the limit d — > °° the FPT probability density for E < is not normalizable since the probability to never pass i = is 
positive. For the E > (and d °°) case it is known Jl^,l3l that the mean FPT distribution density scales as f~' 1+f1 ' for large 
t. This implies that the value of fi may be evaluated by finding the largest converging moment of the mean FPT for d — > °°. 
Equivalently, jj. is given by the smallest moment that diverges 

|i=inf{m:(F")=°o}. (6) 

To obtain an approximation for the full FPT distribution, in the next section we introduce and solve a RT model. The parameters 
of the RT model are set by the known properties of the RF model. 

III. THE RANDOM TILT MODEL 



Recently Oshanin and Redner [23] used an optimal fluctuation method 112711 to study the splitting probability, P, of a random 
walker on a random forcing energy landscape. The splitting probability is defined as the probability to reach site ; = d before 
hitting site i = 0, starting from site i = iq. Within this approach one replaces each specific realization of the random forcing 
energy landscape by its average slope, U, which is specified by the energy difference between the origin and the last, i = d, 
site divided by the total number of the sites (see Fig. |2). The central limit theorem implies that the probability density of the 
average energy tilt is Gaussian: 



iwn : ^- T =S^ ! (7) 
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FIG. 2: On this schematic plot we demonstrate the approach described in t23ll where each specific realization of the disordered is replaced by 
constant energy slope. 



where E, as stated before, is the average energy difference between two subsequent sites 



E = (Ei) = ( In 



Pi-i 



(8) 



and a 2 is the variance of the energy difference between two subsequent sites 



a 2 = (E 2 ) (E,Y 



In 



2 <li 



Pi-l 



ln- 



Pi-i 



(9) 



Using the fact that on a constant energy tilt, U, the splitting probability is given by 1 Il28ll the average splitting probability 



can be well approximated by 12311 



au-ey 



~2^~ 1 - e u '° 



In 



^ \-e Vd 



dU. 



(10) 



Inspired by this approach we introduce a random tilt (RT) model. Within the model the energy landscape of each realization 
is flat while its tilt is a random variable. Namely, defining hopping rates to the right and left by u and v respectively, with 
^ = e~ £ , we take £ to be a random variable with a Gaussian probability density 



Pr(e) 



,(e-(e)) 2 



y/2na 2 /d 



(11) 



with a mean (e) and a variance a 2 . Here d, as before, is the size of the system. As we show below by fixing (e), a £ and 
an overall time we can approximate very well the first passage behavior of a random walker on a RF energy landscape in the 
limit d — > °o and large t . We refer to this as the universal regime. In this paper, however, we are interested in an approximation 
over any time scale when d is finite. In this case the small t and large t behaviors are not universal. As we show a very good 
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approximation can be achieved by taking the value of v as a random variable with its own distribution Pr (v). In Sees. IV Al and 
IV Bl we show that, in order to approximate the FPT probability density of the RF model, the simplest choice one can make is 

Pr(v) = 7Ti5(v- vi) + na8(v — V2) + (1 - %\ - %i)8{v -v^) . (12) 

The values of (e), <7 £ , v\, V2, V3, 7i\ and %2 are set by known properties of the RF model (with no fitting parameters). As stated 
above, the results are summarized at the end of the paper. Note that to approximate only part of the FPT distribution a simpler 
choice of Pr (v) can be made. For example, if one is not interested in the short time behavior one can chose Pr (v) = 8 (v — vo) 
where vq is specified in what follows (see Eq. d33l > below). If one is not interested in the long time behavior but wants to capture 
the short time behavior one can chose Pr (v) = 5 (v — (q)) (see Eq. (|37| > below). 



IV. THE PROPERTIES OF THE RT MODEL 

To approximate the FPT distribution on a RF energy landscape we evaluate the splitting probability of the RT model, Prt, 
the analog of the exponent ji of the RT model, denoted by [Irt and the average FPT, denoted by (t) RT . In addition we also 
calculate the Laplace transform of the FPT probability density of the RT model, Frj (s). 

The splitting probability of the RT model can be easily obtained, similar to Eq. ([Tol l, and is given by 



Prt I - T^f T^fde. (13) 
2*% 



1 - e £, o 

T 



Note that Eq. (fT~3T > is exact for the RT model. 

First we demonstrate that an analog of /x, denoted by Urt, exists for the RT model when rf->°o. This shows that indeed the 
model can reproduce a qualitative first passage behavior similar to a random walker on a RF energy landscape. To do this we 
calculate the scaling behavior of the moments of the average FPT, f £ . For a given realization of e (the equivalent of a random 
realization of the RF energy landscape) this scales as 

1 (n) d - 1 \ e- £d £ < 
h = ~±*J— r ~{ . , • (14) 

v v- 1 £>0 



Recalling that the probability density of e is Gaussian 



v 1-e" 



Pr(e)~e 2ct e , (15) 

and in analogy with the RF model (see Eq. ©), [Irt is given by 

PLrt = inf {m : = <*>} , (16) 

where now the angular brackets denote an average over £ values. In the d — > °° limit, when (f7' n ) = °°, the average of the m'th 
moment is controlled by the contributions where e < so that 

<£)~A ^) . (17) 
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We, therefore, obtains 



2(e) 



(18) 



Thus, as stated above, an analog of ji exists within the RT model after averaging over realizations of disorder. It, of course, 
does not exist for each realization of the model. Note that this expression is identical to that obtained for a full random forcing 
model where the random force is drawn from a Gaussian distribution with a mean (e) and variance Jl^j . Furthermore, Eq. 
(TT~4t implies that the scaling behavior of the typical FPT is given by 



exp (\nt) RT ~ e 



-{e)d 



(19) 



increasing exponentially with d for (e) < and remaining constant when (e) > 0. 

Next we consider the mean FPT of the RT model (averaged over realizations of the tilt of the landscape) to the origin from 
site iq. For a given energy tilt, £, and a given left hopping rate, v, the thermal averaged FPT may be obtained using Eq. ® with 
the substitutions ( |: / = v anc ^ ( e ~ E ') = e& - Averaging over e and v one obtains 



(J) 



RT 



(20) 



where 



/ = 



j-e(d+l) _ g -e(</-io+l) 
(e- £ -l) 2 



'0 



-1 



2** 



-de 



(21) 



and 



-)= fipr(v)dv. 
v Jo V 



(22) 



The Laplace transform of the FPT probability density can also be obtained exactly and is given by 



Frt (s) 



Vr{e)Vx(v)<S>(-,e\dedv 



(23) 



where <t> (^,e) is the FPT probability density on a flat energy landscape with a tilt e and a left hopping rate v. Using standard 
FPT results B28I1 one has 



<F 



2'» 



(24) 



with 



(25) 



Next we use these results to approximate the FPT distribution in the RF model. 



V. APPROXIMATING THE FPT DISTRIBUTION OF THE RF MODEL 



To use the results to approximate the FPT probability density of the RF model here we set the parameters (e) and a £ and 
the probability distribution Pr (v), such that Frt (s) yields a good approximation to F (s). Namely, we work with the Laplace 
transform. The matching is done so the distributions agree in different ranges of s and can be done by matching between known 
quantities of both the RT and the RF models. We first determine the parameters (e) and a e . 

A. Determining (e) and a £ 

Here we find expressions for (e) and a £ by matching the scaling properties of the FPT probability densities. The value of 
(e) may be found by matching the scaling behavior of the typical FPT of the two models in the large d limit. We showed above 
(see Eq. Q) that in the RF model the typical FPT grows as e~ Ed for E < and remains constant while E > (see Eq. (fT4b). 
The same is true (see Eq. (fT9b ) for the RT model when E is replaced by (e). Thus to match the scaling behavior of the typical 
FPT in both models we set 

(e)=E. (26) 

The scaling properties in the intermediate s regime for a finite value of d can be expected to behave identically to the small 
s behavior of the d — > °° limit. Since the scaling behavior is different for negative and positive average tilts we separate our 
discussion to two cases. 

1. The case E < 

In this case the splitting probability, P, remains finite in the d — > °° limit. Thus, the small s behavior of the FPT probability 
density in the d — > °° limit is 1 — P + O (s). To match the intermediate s behavior of FPT probability density of RF and RT 
models we match the splitting probabilities (i.e. the leading, O , order of the FPT probability densities) of these two models. 

Comparing ( TTcTb with ( fTST l and using Eq. (l26T i yields 

CT £ = C7. (27) 
Note that both the average tilt and the variance in the RF and RT models are equal in this regime. 

2. The case E > 

As studied in Sec. Q] in this case the mean FPT distribution for d — > °° in the small s region scales as 1 — O (s^) for /i < 1 
(or, equivalently in time as t~^ l+ ^). Therefore, we first match /Xrt of the RT model with fx of the RF model. Namely, we set 
Urt = Using Eq. (fT8T i this yields 
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Solving Eqs. (f26t and (l28l l we get 




Note when the distribution of £, is Gaussian equations ( |27| > and ( |29] l become identical. Namely, both the average tile and the 
variance of the RF and RT models are equivalent. When the distribution is not Gaussian the average tilt is the same in both 
models but the variances can be different. In this regime it is important to match the precise value of fi rather than the splitting 
probability as in the E < case. 

B. Determination of Pr (v) 

In this section we choose the distribution Pr(v) by matching the FPT probability density of the RT model to that of the 
original RF system in the small and large s limits. We show that for each limit one may choose v to be a constant. However, this 
constant depends on the limit. Thus, to match the full s behavior v cannot be a constant and, as we show below, the simplest 
choice for Pr (v) is 

Pr(v) = %5(v-vi) + %5(v-V2) + (l - %\ - wj) 8 (v- v 3 ) (30) 

with Vi t 2,3 are set according to relevant time scales in the problem and 7Ti,2 are determined by the matching procedure, described 
below. 



1. The small s approximation 
The s — > behavior of the FPT probability densities of the RF and the RT models is 

F(s^0) = 1 - (t)s (31) 



and 

F RT (s^0) = l-{t) RT s, (32) 

respectively. Therefore, to match F (s — > 0) and Frt (s — > 0) one has to match the mean FPTs of both models. Demanding 
(J) RT = (J) and using Eq. d20l one has 

OK 

where / is defined in Eq. (l2Tb and (J) is given by Eq. ©. 

Below we use this fact to match the behavior of the FPT distribution of the RF model over the whole s range. Note however, 
that if one is not interested in the short time behavior the procedure is simplified. One can then choose Pr (v) = 8 (v — Vi) where 
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■j|y. Then, using Eq. (l23l . the Laplace transform of the FPT probability density is given by 



2ai i s 



F RT {s)= _ 4> — ,e )de. (34) 



2. 77ie /arge s approximation 

The large 5 limit corresponds to the small t limit. This regime of the FPT distribution is controlled by walks which hop only 
to the left before reaching the origin. Therefore, in the s — > °° limit the behavior of the FPT probability densities of the RF and 
the RT models is 

F(*->oo) = < 9 >*> (35) 

and 

F RT (s^oo) = (v i '>), (36) 

respectively. Thus, to match the large s behavior of Frt (s) and F (s) we demand 

(v io ) = (q} i °. (37) 

Next we match the behaviors in the RT and RF models over the whole range of s. However, we comment that if one is not 
interested in the approximate FPT distribution in the long time limit one can choose Pr(v) = 8 (v — V2) where Vi = (q). Then, 
using Eq. d23l . the Laplace transform of the FPT probability density is given by 

F RT {s)= _ <P[-,e)de. (38) 



2*$ 



V 2 



3. An approximation for the whole s range 

To match both the small and large s limit one should supply Pr(v) such that Eqs. d33l and d37b are satisfied. However, these 
equations do not determine Pr (v) uniquely and there is a lot of freedom. As stated above, the simplest choice that satisfies these 
equations is 2 

Pr(v) = n\8(v- vi) + n 2 8(v- v 2 ) + (1 - n\ - n 2 ) 8 (v - v 3 ) . (39) 

However, since one has to satisfy only two equations ( d33l and d37] >) there is still a lot of freedom with five free parameters. 
Therefore, first we set Vi 2 3 such that they represent relevant time scales in the problem. Eq. (1331 suggests 

i-f. (40) 

VI / 



: Other, simple choices lead to nonlinear equations. For example, using Pr(v) = |5(v — Vi) + |5 (v — V2), Eqs. {33) and )37t lead to ~ + ^ = 2-*f- and 

v'°+v^=2{ 9 )'». 
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while Eq. (|37| > suggests 

V 2 = (q). (41) 

The third time scale is chosen to represent the average time of a hop to the left 

\>3 = (-\ ■ (42) 



,<7 

Given these three quantities, vj 2.3, Eqs. (f33T > and d37l i give 



n lV \° + 7C 2 v'° + (1 - m - % 2 ) v 3 ° = (q) h (43) 



and 



Vl V2 V 3 / 



(44) 



The solution for Tt\ 2 is 



V2V3 I V 2 — V 3 j + Vj (V2 — V3J + V1 I v 3 — v 2 I 



and 



„ _ v 2 (r) v ' v H V ' - y 3j+i^(^-^) + 7 y(v3 -Vj , 
^2 — T 7 : 7T : ~ ; 7 : ~ ; : , , \ • 



V 2 V3 I V 2 - V3 I + Vj (V 2 - V3J + V'! I v 3 — v 2 I 



Using these quantities, the probability density (f39b for v and Eq. (1231 . the Laplace transform of the FPT probability density is 
given by 



Frt (s) = 



,t,4>( ^-,ej +% 2 ®{^e\ +(l-7t2-7c 3 )<$( -^,e 



de. (47) 



In the next section we compare the obtained result to the numerical simulation. This expression with the above chosen values 
of parameters is the main result of our paper. Note that while quite a few parameters have to be set there are no free fitting 
parameters. As we show the agreement is very good. 

C. FPT probability density and comparison to numerical result 

In this section we compare between Eq. ( f4Tb and numerical simulations. Namely we compare the Laplace transform of the 
FPT probability, density, Frt (•?), and the survival probability 

F s {t) = \- f F(t')dt' (48) 



of the RF and RT models. We checked the results for RF models with both Bernoulli and Gaussian disorders. 
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FIG. 3: In this graph a comparison between the analytic approximation d47t and the numerical results for a Bernoulli disorder is shown in 
Laplace (top) and time (bottom) spaces. The parameters for these plots are: r = 0.5, £\ = —0.4, £2 = 0.4 and d = 100, such that jj. = 0. 
Circles represent the choice ;'o = 1 while squares represent the choice i'o = 50 cases on both, left and right, figures. The red lines are analytic 
approximations based on Eq. d47l >. 

For the RF system with Bernoulli disorder the energy difference on each site, is drawn from the distribution: 

r Ej = £1 

Pr (£;) = < ■ ( 49 ) 

1 - r Ei = e 2 

On Figs. !3l4l5l one may see the comparison between Eq. d47T > and the numerical calculation in different parameter ranges. As 
can be seen the agreement is very satisfying. 

For a RF model with a Gaussian disorder the energy difference between sites, is drawn from a normal distribution with a 
mean E and a variance a 2 : 

Pr(£,) = ^=^. (50) 

On Figs. I6l7l8l we show a comparison between Eq. d47l > and numerical calculation in different parameter ranges. Again, the 
results of the approximation are very good. 



VI. SUMMARY 



In this paper we presented a random tilt model and solved it analytically. We showed that the model can be used to obtain an 
approximate expression for the FPT distribution of a random walker on a random forcing energy landscape. To do this several 
parameters of the RT model have to be set as a function of the RF model parameters. As we showed, this can be done with no 
free fitting parameters. For convenience we summarize the results below. 
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FIG. 4: In this graph a comparison between the analytic approximation d47t and the numerical results for a Bernoulli disorder is shown in 
Laplace (top) and time (bottom) spaces. The parameters for these plots are: r = 0.49, £j = —0.1, £2 = 0.1 and d = 550, such that fj, = 0.4. 
Circles represent the choice io = 1 while squares represent the choice z'o = 100 cases on both, left and right, figures. The red lines are analytic 
approximations based on Eq. d47l >. 




FIG. 5: In this graph a comparison between the analytic approximation d47t and the numerical results for a Bernoulli disorder is shown in 
Laplace (top) and time (bottom) spaces. The parameters for these plots are: r = 0.58, £\ = —0.3, £2 = 0.3 and d = 150, such that fl = —1.08. 
Circles represent the choice ;'o = 1 while squares represent the choice (q = 50 cases on both, left and right, figures. The red lines are analytic 
approximations based on Eq. d47| >. 
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FIG. 6: In this graph a comparison between the analytic approximation d47t and the numerical results for a Gaussian disorder is shown in 
Laplace (top) and time (bottom) spaces. The parameters for these plots are: E = 0, O" = 0.3 and d = 150, such that fl = 0. Circles represent 
the choice i'o = 1 while squares represent the choice iq = 50 cases on both, left and right, figures. The red lines are analytic approximations 
based on Eq. ( I47t . 




FIG. 7: In this graph a comparison between the analytic approximation d47t and the numerical results for a Gaussian disorder is shown in 
Laplace (top) and time (bottom) spaces. The parameters for these plots are: E = 0.01, o" = 0.3 and d = 250, such that fj. = 0.2222 Circles 
represent the choice i'o = 1 while squares represent the choice i'o = 50 cases on both, left and right, figures. The red lines are analytic 
approximations based on Eq. d47| >. 
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FIG. 8: In this graph a comparison between the analytic approximation d47t and the numerical results for a Gaussian disorder is shown in 
Laplace (top) and time (bottom) spaces. The parameters for these plots are: E = —0.01, (7 = 0.2 and d = 350, such that fl = —0.5 Circles 
represent the choice io = 1 while squares represent the choice i'q = 50 cases on both, left and right, figures. The red lines are analytic 
approximations based on Eq. ( 1471 ). 



The approximation for the random walker's Laplace transformed FPT probability density from site io > to the origin on a 
random forcing energy landscape with i.i.d. random variables {/?,} and } is: 



F(s) = 



2*4 



7Ti4>( — ,£ 



7Z 2 <t>( — ,£ ) + (1-7T2-7T3)* ( — ,£ 



de (Eq. (|47]i), 



where 



4> 



i •{) 

A 2 



2'o 



MP) 



d-i'o-l 



MM 



d-1 



,e -2 



(f + l)^(f, £ )-2 



(Eq. (El), 



^ 1 ' 2 ^ 1 + + ~) ' ± ^l+^(J-l) +e -(J + l) 2 (Eq. 



£3), 



W {{e- E -)-lf \ qi /{e-*i)-l 



(Eq. ©), 



15 



/ = 



,-e(d+l) _ e -e(d-io+i) 
(e- £ -l) 2 
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, (e-(e)) 2 



2k4 



-de (Eq. (ED), 



V, = ^=y (Eq. (S3), 



v 2 = ( 9( > (Eq. (ED), 



-l 



v 3 = ( - ) (Eq. (03), 



(F)V 2 V3(v|-V^)+^(v 2 -V 3 ) + ^( 1 



,'o+l ,,'o+l 



/ f '0 '0 1 i '0+1 C \ I ( in+1 'o+l \ 

v 2 v 3 - v 3 u J + v{' (V2 - v 3 ) + vi (^v 3 u - v 2 u j 



(Eq. 65J), 



^2 : 



V2 



(F) v 1 v 3 (v'r^)+^(v 1 -v 3 ) + ^(v: 



,'0+1 ,,'o+l 



/ ( 'o k>\ i 'o+l z' \ i f i'o+i io+l \ 

V 2 V'3 I V 2 U - V 3 U I + Vj" (v 2 - v 3 ) + vi I v 3 u - v 2 u ] 



(Eq. ©), 



(e) = £ (Eq. and 



a £ <0 



£ >0 



(Eqs. (|22ll and ((29j). 



If one is not interested in the large or in the small s behavior one may use the much simpler Eqs. d34i > or (l38l >. respectively, 
instead of Eq. < f4Vl >. 

Comparing this approximation with the numerically calculated FPT distribution of the RF model we showed that the first 
may serve as a good approximation to the second. Finally, similar methods can be used to approximate other, say absorbing, 
boundary conditions at i = d + 1 . 
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